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Abstract. We present a fully relativistic and gauge invariant framework for calculating the cross 
sections of dilepton production in nucleon-nucleon (NN) collisions which is based on the meson- 
exchange approximation for the NN scattering amplitudes. Predictions of our model are compared 
with those of other covariant models that have been used to describe this reaction. We discuss the 
comparison of our calculations with the old DLS and the recent HADES data. 
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1. Introduction 

Experiments with High Acceptance Di-electron Spectrometer (HADES) are aimed at 
searching for medium modifications of hadrons at moderately high temperatures (T < 100 
MeV) and baryonic densities up to 3 times the normal nuclear matter density. Due to neg- 
ligible final state interactions with surrounding medium, dileptons (e+e^) provide a very 
clean and powerful probe for this purpose. 

A recurring feature of the dilepton spectra measured in nucleus-nucleus (AA) colli- 
sions has been the enhancement (above known sources) in the invariant mass distribution 
of the cross sections in the intermediate region of dilepton invariant masses (M). This 
has been the case for experiments performed for bombarding energies ranging from as 
low as 1.0 GeV/nucleon (DLS and HADES data [1,2]), through the SPS energies (40 - 
158 GeV/nucleon) [3-5] to the energies employed by the PHENIX collaboration at RHIC 
(which correspond to invariant mass of 200 GeV/nucleon) [6]. The enhancement seen 
at the SPS energies are understood in terms of the modification of the p meson spectral 
function in the hadronic medium [7]. 

The large dileptons yields observed in the DLS experiment (even for the light collision 
systems ^^C + ^^C) in the invariant mass range of 0.2 - 0.5 GeV, are yet to be explained 
satisfactorily [8-12]. Independent transport model calculations (TMC) have been unable 
to describe these data fully even after including contributions from (i) the decay of p and 
u! mesons which are produced directly from the nucleon-nucleon (NN) and pion-nucleon 
scattering in the early reaction phase [13], (ii) the in-medium p spectral functions [14], (iii) 
the dropping p mass with coiTesponding modification in the resonance properties [8], (iv) 
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an alternative scenario of the in-medium effect - a possible decoherence between the inter- 
mediate meson states in the vector resonance decay [11]. This led to term this discrepancy 
as "DLS-puzzle" [8,13] which persists even now. The new dilepton production data of 
the HADES Collaboration agree remarkably well with the corresponding DLS data [15]. 
Therefore, there is no longer any question against the validity of the DLS data and the 
failures to explain them by various transport models have to do with problems in the theo- 
retical calculations. 

Cross sections for dilepton production in elementary nucleon-nucleon (NN) colUsions 
are important inputs in the TMC. For the NN bremsstrahlung processes most TMC use the 
cross sections predicted by a semiclassical soft photon approximation (SPA) model [16,17]. 
An important point made in a recent TMC [18], is that if one scales up the SPA cross sec- 
tions for the pp and pn systems by factors of 2-4, the observed dilepton yields of both DLS 
and HADES experiments can be reproduced. These enhancement factors are motivated by 
the larger elementary TV TV dilepton cross sections predicted by the calculations reported 
in Refs. [19,20]. However, these large cross sections are not in agreement with those re- 
ported in several previous studies [21-23]. They also do not agree with results of a more 
recent study of the dilepton production in NN collisions [24] . Therefore, the puzzle seems 
to have reduced to a proper understanding of the dilepton production in elementary A''A'^ 
collisions. 

In the next section we present a brief comparison of various models used to calculate the 
dilepton yields in elementary NN collisions. It is important to understand the differences 
seen in the predictions of various models in order to have a proper theoretical description 
of the new HADES data on the dilepton production not only in elementary pp and pn 
collisions but also in the AA collisions. 



2. Comparison of various models of dilepton production in elementary NN Collisions 

Calculations reported in Refs. [19,20] (to be referred as KK), and [23] (to be referred as 
SMI) use the same basic model to calculate the dilepton production in NN collisions. 
Both KK and SMI account for the initial interaction between two incoming nucleons by 
an effective Lagrangian which is based on the exchange of the tt, p, co and a mesons. The 
couphng constants and the form factors at the nucleon-nucleon-meson vertices are the same 
in both studies. These parameters were determined in Ref. [21] by directly fitting the T 
matrices of the NN scattering in the relevant energy region, and they have been used in the 
successful descriptions of the NN NNn [25], pp pAK+, pp pTPK+ [26,27] 
and NNrj [28] reactions. 

The major difference between calculations of Refs. [19] and [23] lies in the method of 
implementing the gauge invariance of the NN bremsstrahlung amplitudes. To investigate 
this issue, the dilepton production cross sections within the SMI model have been recalcu- 
lated recently [24] (to be referred as SM2) by using a pseudoscalar (PS) coupling 

CnNt, = igNNTv'^NlbT ■ *ir*Ar, (1) 

for the nucleon-nucleon-pion (NNtt) vertex while the KK model employs a pseudovector 
(PV) coupling 
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Figure 1. The invariant mass distribution of the A'^A'^ bremsstrahlung contributions to 
the dilepton spectra in pp (left) and pn (right) collisions at the beam energies of 1.04 
GeV and 2.09 GeV. Results obtained within our model are shown by solid lines while 
those of Refs. [22] and [19] by dotted and dashed lines, respectively. 

for this vertex. With a PV vertex an extra term (known as contact or seagull term) ap- 
pears in the model which comes from the electromagnetic coupling (by the substitution 
— iemA^, where m is +1, 0, -1 for positive, neutral and negative pions) in the 
NNtt Lagrangian. On the other hand, with a PS vertex, the theory is free from this term. 
For point-Uke particles both prescriptions lead to gauge invariant NN bremsstrahlung am- 
plitudes. However, nucleons have finite sizes - to account for this fact form factors are 
introduced at the corresponding vertices. Form factors are also required to quench the un- 
physical contributions from higher momenta. In the presence of form factors the gauge 
invariance may be violated. As is shown in Ref. [21], in case of the PS coupling it is 
straight forward to introduce the form factors while retaining the gauge invariance by fol- 
lowing the method of Ref. [29]. However, with contact terms present in the theory in the 
PV case, the introduction of form factors without violating the gauge invariance, is a more 
delicate issue. There are several prescriptions for this and various procedures may lead to 
quite different results (see, e.g., [30]. Therefore, the use of the PS coupling for the NNtt 
vertex makes the dominant pion exchange contributions free from such ambiguities. 

In Ref. [22] (to be referred as dJM model), instead of the one-boson exchange picture 
of SMI, SM2 and KK calculations, the nucleon-nucleon interaction is included via T- 
matrices that are based on the Paris potential. This model provides a reliable description of 
the off-shell behavior of the effective NN interaction. In addition, the rescattering terms 
are included in both nucleon and A intermediate states - these terms are not present in both 
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Figure 2. The invariant mass distribution of ttie A isobar contribution to the dilepton 
spectra in pp (upper panel) and pn (lower panel) collisions the beam energy of 1.04 
GeV. The results of our model are shown by full lines while those of Ref. [19] by 
dashed line. 

SM (1 & 2) and KK calculations. Furthermore, the degeneracy of the PS and PV choices 
for the pion coupling is resolved to a large extent in this model; one has to simply use 
different parameters to obtain the same fit. However, the price paid here is that the nucleon 
current is not gauge invariant. These authors rectify this problem in an ad-hoc manner. 

In the next section we present a comparison of the predictions of various models for the 
invariant mass distributions of the dUeptons produced in NN collisions. 

3. Comparison of cross sections calculated within various model 

In Fig. 1, we show the invariant mass distribution of the pp and pn bremsstrahlung con- 
tributions to the dilepton spectra at beam energies of 1 .04 and 2.09 GeV as calculated in 
SM2 model (sohd lines). Also shown here are the results for this reaction as reported in 
Refs. [22,19]. We first discuss the results for the pp reaction shown in the left panel of 
this figure. At the outset we remark that the SM2 cross sections are very similar to those 
reported in Refs. [21] and [23]. However, some differences are seen between the cross sec- 
tions of SM2 and dJM models and also between them and the KK calculations. Whereas 
the dJM cross sections are larger than those of SM2 in the mass region M < 0.25 GeV, 
they are almost identical to the latter at 2.09 GeV in this mass region. At both the beam 
energies, dJM results are smaller than SM2 ones for M > 0.25 GeV. In contrast to this, 
the KK cross sections are larger than those of SM2 everywhere for M < 0.6 GeV. An 



4 



Pramana - J. Phys., Vol. xx, No. x, xxx xxxx 



CM 

o 

> 

0) 

O 



■a 



10 



1 


1 ' 


1 ' 1 ' 1 ' 1 ' 1 

■ ■ ■ ■ % decay I 


r 










— — Kaptari & Kaempfer : 


- 

r 




, ^ 1 .Ut vac V -a 




//^ 






if- 


JS^ PP = 




11 \ 










J 


i 1 , 


1,1.1. i7- 



0.1 0.2 0.3 0.4 0.5 
M (GeV/c^) 



Figure 3. The calculated dilepton invariant mass distribution for the pp collision at the 
beam energy of 1.04 GeV in comparison to the DLS data. The contribution of the tt'^ 
Dalitz decay is also shown here which is the same as that in Ref. [23]. 



important point to note is that there is no overall multiplicative factor that differentiates the 
results of various models. 

It is rather surprising that despite using the same diagrams, input parameters and gauge 
invariance restoration procedure, the SM2 pp bremsstrahlung cross sections are lower than 
those of Ref. [19]. Of course, in Ref. [19] a pseudovector NNtt vertex has been used as 
compared to the pseudoscalar one employed in SM2. In this context, it is worthwhile to 
note that for the real photon production, the covariant model calculations do not depend on 
the choice of the NNtt coupling (PS or PV) as is shown in Ref. [31]. In case of dileptons, 
different results can arise for two couplings from the magnetic part of the NNj vertex. In 
fact, in Ref. [32] it is shown that pp dilepton bremsstrahlung contributions obtained with 
the PV NNtt coupling are actually smaller than those calculated with the PS one at the 
beam energy of 2 GeV. The calculations presented in Ref. [23] also support this to some 
extent. 

In the right panel of Fig. 1, we compare the SM2 results with those of the KK model 
for the pn collisions. In this case, the situation is even more contrasting - the SM2 cross 
sections are a factor of 3-4 lower than the KK ones. The results of the dJM model are not 
available for this reaction. It is highly desirable that this discrepancy between the results 
of the two models for the pn case is properly understood. 

In Fig. 2, we show a comparison of SM2 and KK results for the invariant mass distribu- 
tion of the A isobar contribution to the dilepton production in pp and pn collisions at the 
beam energy of 1 .04 GeV. We note that here too the KK cross sections are larger than SM2 
by factors of ^ 2 at smaller values of M even though the two models have used the same 
ingredients and input parameters for this part and there is no ambiguity related to gauge 
invariance as the resonance vertex is gauge invariant by its very construction. 

In Fig. 3, we compare the SM2 and KK total cross sections for the dilepton production 
in pp colhsions with the DLS data at the beam energy of 1 .04 GeV. The cross sections 
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Figure 4. Contributions of n (dotted line) and p (dashed line) exchange diagrams to 
the invariant mass distribution of the dileptons produced in proton-neutron collisions at 
the beam energy of 1.25 GeV. Their coherent sum is shown by the solid line. 



calculated within SM2 model are folded with appropriate experimental filter and final mass 
resolution. The folded KK cross sections have been obtained by assuming that the folding 
procedure does not affect the ratios of the unfolded cross sections in the two cases. In 
this figure we have also shown cross sections for the 7r° Dalitz decay (7r° 76+ e~) 
which are the same as those in Ref. [23]. It is seen that KK cross sections overestimate 
the DLS data for M < 0.3 GeV where statistical errors are smaller. The data have larger 
error bars for M > 0.3 GeV. In this context, the HADES data on the elementary dilepton 
production reactions are useful because of their low statistical error. A comparison of the 
HADES dilepton mass spectra measured in the pp and quasi-free np reactions [33] and the 
predictions of the SM2 and KK models, has been reported in Ref. [34]. There one notices 
that for the case of pp reaction, once again the predictions of the SM2 model are in pretty 
good agreement with the data while the KK model overestimates them. 

In Fig. 4, we examine the role of various meson exchange processes to the dilepton 
production in the pn collisions at 1.25 GeV beam energy. In this figure, the dotted and 
dashed lines show the individual contributions of the n- and p-exchange graphs to the 
dilepton invariant mass distribution. Contributions of other meson exchange processes are 
relatively smaller. The coherent sum of the two terms is shown by the solid line. It is 
clear that the 7r-exchange process dominates the cross sections and more significantly, it 
interferes destructively with the p-exchange term as the coherent sum of the two diagrams 
is even smaller than the pion-only exchange term. This impUes that a larger p-exchange 
term [20] will lead to even smaller total cross sections. 
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Figure 5. (Left panel) the dilepton invariant mass distribution in np (dashed line) 
and quasi-free np interaction (solid line) at 1.25 GeV/nucleon. In the later case the 
center mass energy is smeared by including the neutron momentum distribution in the 
deuteron using the Argonne VI 8 potential. (Right panel) the dilepton invariant mass 
distribution in the quasi-free np collisions at 1.25 GeV/nucleon calculated within the 
model of Ref. [24] (dashed line). Also shown are the contributions of ri meson and 
A(1600) Dalitz decay processes by dashed-dotted and dotted lines, respectively. Their 
incoherent sum is shown by the solid line. 



3.1 Calculations at the HADES energies 

In the experiments carried out by the HADES Collaboration, proton and deuteron beams 
with kinetic energies of 1.25 GeV/nucleon were incident on a proton target. In this mea- 
surement the quasi-free np reactions were selected by detection of fast spectator protons 
from the deuterium breakup in the forward direction. 

As has already been stated above, the HADES data on the dilepton invariant mass distri- 
bution (DIMD) in the pp reaction were described well by the SM2 model [33,34] while the 
KK model overestimated them. However, for the quasi-free np collisions, both models fail 
to reproduce the experimental cross sections. The shape of the DIMD changes dramatically 
when going from pp to pn interactions. One notices that (1) in the dilepton invariant mass 
region of 0. 15 to 0.35 GeV, the yield is about 9 times higher in the np case as compared to 
that of the pp case, (2) in the quasi-free np case, the tail of the DIM extends to much larger 
values of M, (3) the ratio of quasi-free np to pp cross sections reaches almost a value of 
100 at M = 0.5 GeV. 

To understand points (l)-(3), we introduce the following 3 additional features in the SM2 
model for the np case, (i) The available energy in the center of mass is smeared to include 
the neutron momentum distribution in the deuteron using the Argonne V18 potential [35]. 
The consequence of this procedure is that the dp reaction results is a smeared np reaction 
with center of mass energies in excess of the threshold for the 77 meson production (see, 
e.g., Ref. [36]). (ii) Hence, the contribution from the 77 Dalitz decay is taken into account. 
In these calculations the total cross sections for the n+p^n+p^-q reaction is taken 
from the Ref. [28] where a good description is achieved of the available corresponding 
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experimental data, (iii) The dilepton yields from the Dalitz decay of higher mass delta 
resonances are included. Our preliminary results are shown in Fig. 5. It was already noted 
in Ref. |34| that the predictions of the SM2 model are in agreement with the np data for 
M < 0.3 GeV. In the left panel of Fig. 5, we note that the Fermi folding procedure leads 
to an extended tail in the invariant mass distribution. In the right panel, we show the total 
cross sections obtained by Fermi folding the cross sections of the SM2 model together with 
contributions of the rj meson and A(1600) isobar Dahtz decay processes. Their incoherent 
sum is shown by the solid line. It is seen that the inclusion of contributions (ii) and (iii) 
leads to only a marginal enhancement in the summed cross sections over those of the SM2 
model for M values around 0.50 GeV. Therefore, a proper understanding of point (3) is 
stiU lacking at this stage. 

4. Summary and Conclusions 

In summary, we have shown that the dilepton yields from the NN bremsstrahlung pro- 
cess calculated within the SM (1 & 2) models [23,24] are lower than those of the KK 
model [19,20] by factors of 2-4 for both pp and pn collisions at 1 - 2 GeV incident ener- 
gies. The recent HADES measurements have confirmed that the cross sections calculated 
within the SM2 model are in good agreement with the experimental data for pp reaction 
while those of the KK model overestimate them. This was also the case for the DLS data 
as shown in Fig. 3. This impUes that care must be exercised in using the larger NN 
bremsstrahlung cross sections in the transport model calculations in order to explain the 
DLS and HADES data obtained in the nucleus-nucleus colhsions. 

With HADES collaboration confirming the old DLS data, the resolution of the "DLS" 
puzzle has indeed now shifted to a better understanding of the theoretical models. In this 
context, the focus should be to achieve a proper explanation of the HADES data for dilep- 
ton production in quasi-free np interaction. We showed that smearing of the center of mass 
energy to include the neutron momentum distribution in deuteron, does lead to extended 
tails in the dilepton invariant mass distributions. However, the absolute magnitudes of the 
cross sections in the region of dilepton invariant masses around 0.5 GeV are still not ex- 
plained even after taking into account the contributions from the Dalitz decays of rj meson 
and higher mass resonances. Work is in progress to resolve this new puzzle. 
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